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ABSTRACT 


Mean flow and disturbance measurements have been made on a flat 
plate model, one meter long in the CWRU low speed wind tunnel using 
hot-wire anemometry. Tests were conducted on a smooth plate and with 
two different grades of sandpaper covering; first with a sandpaper 
covering giving Re^ values centered around 10 and then with a coarser 
sandpaper having Re^ values of about 150. 

The mean flow measurements indicated Blasius profiles for the 
smooth plate and with the finer roughness. With the coarser rough¬ 
ness, the profile was best approximated by a Blasius profile at the 
early stations but then progressed toward a turbulent profile. 

The finer surface roughness (Rc^ 1 10) had only a very slight 

effect on the u' distributions relative to those for a smooth 
rms 

plate. This showed up primarily through examination of the distri¬ 
bution of local, turbulence intensity, u' /u. which proved to be a 

rms 

much more sensitive indicator than dimensionless turbulent energy, 

2 2 

u' /!' . for the coarser roughness (Re, ~ 150), local u' intensi¬ 
ons c k rms 

ties in the boundary layer of the order of 5% were observed before 
the mean profile departed from a Blasius shape. The progression 
toward turbulent mean velocity profiles was accompanied by large 
increases in turbulent intensities over the whole boundary layer but 
especially near the wall. A study of disturbance spectra for this 


i i 






<';isc shows largest amplitudes and amplifications at frequencies well 
below the Blasius neutral curve. The maximum amplification rates 
are about three times those for Tollmien-Schlichtinfi instability 
over a smooth plate. Transition occurs at Re ' 350,000 which is 

A 

much below that for a smooth surface. 
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CHAPTER I 


INTRODUCTION 


On a smooth surface, transition of the boundary layer presumably 
occurs as the result of the amplification of the disturbances ini¬ 
tially present in the flow. In this model of the transition process, 
the disturbances amplified are those whose frequencies lie within 
the neutral curve for Tollmien-Schlichting waves and are dependent 
on the Reynolds number of the boundary layer. If transition on a 

smooth plate occurs at distance X downstream of the leading edge 

c o 

(with corresponding Reynolds number. Re = UX /v), then for 

C o 

roughness, transition occurs at some location less than , 

o 

corresponding to a transition Reynolds number, Re^ = UX^/v,which is 

smaller than Re (Fig. 1). Roughness enhances friction and drag, 
o 

and also increases the momentum thickness of the boundary layer. 

A minimum critical roughness is sometimes defined as one for which 

Re is approximately 0.95 Re [1]. 

o 

The treatment of roughness is divided into consideration of 
single roughness elements, both two-dimensional and three-dimensional, 
and of distributed roughness which is considered to be a collection 
of three-dimensional protuberances with statistical distribution of 
height, k, and distributed over the plane of the surface. 

There are many experimental results for single roughness ele¬ 
ments and there is even a basic understanding of the destabilizing 


1 






2 


mechanism. On the other hand, much less has been determined about 
distributed roughness and there is no understanding of its destabi¬ 
lizing mechanism. 

For single two-dimensional roughness elements, Dryden jl] 

presents a curve of Re /Re vs. -rr as the best correlation of the 

t t <S* 

o 

body of data to that time. Smith and Clutter [2] argue that this 
is too simplistic a description since there are also effects of the 
shape of the element and the length of laminar flow ahead of the 
element. They argue that a more fundamental parameter for effects 

v 

in the vicinity of a roughness element is Re^ = —. They found 

that for a given shape of roughness. Re, . . , was found to be 

k,critical 

relatively constant and independent of roughness height and location, 
pressure distribution and free-stream turbulence level. There was, 
however, a major effect of roughness shape. 

A most convincing argument on the mechanism by which a two- 
dimensional roughness element destabilizes a laminar boundary layer 
is by Klebanoff and Tidstrom [3]. They show, based on their experi¬ 
ments, that the basic mechanism by which a roughness element induces 
earlier transition is through the destabilizing influence of the 
wake-like recovery zone of the element (Fig. 2). Disturbances in 
the boundary layer are greatly amplified in the recovery zone. The 
growth rates are many times those of the unaffected Blasius layer. 

The fact that the local profiles in the recovery zone depend prima¬ 
rily on the nature of the roughness element establishes Re^ as the 
primary parameter of the phenomenon, and the dependence on the pro- 









perti.es of the parent flow is very weak. 

For three-dimensional roughness elements, the critical values 
of Re^ are higher than for two-dimensional roughness but the desta¬ 
bilizing mechanism may be similar in that disturbances can be great¬ 
ly amplified by the wake-like unstable profiles downstream of the 
elements. 

In view of the results for single roughness elements and the 
known sensitivity of boundary layer stability to small profile 
changes, the analytic investigations to date of distributed roughness 
are all based on the premise that distortion of the mean flow by the 
roughness enhances instability and leads to earlier transition. 

Singh and Lunfley (4 ] made calculations to find the roughness 
induced mean velocity profile which they expressed in terms of the 
spectral density of the roughness and an influence function. The 
resulting profile had an inflection point deep in the inner viscous 
region and in their opinion does not contribute to instability. 
Outside the wall region, the profile they obtained is more stable 
since the curvature in the vicinity of the critical layer has been 
increased. Since experimentally, roughness leads to earlier tran¬ 
sition, Singh and Lumley argue that even if the profile is more 
stable, the roughness also introduces disturbances with energy at 
wavenumbers to which the profile is unstable. 

Lessen and Gangwani [5] model the problem by introducing a wavy 
wall to represent one component of the roughness spectrum. The 
equation for this roughness induced disturbance is the Orr-Sommerfeld 
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equation with phase velocity equal to zero, but with a non-homoge- 
neous boundary condition at y = 0. They retained the Reynolds 
stress term in their mean flow equation and so calculated the dis¬ 
torted mean flow. Once this distorted profile was found, the stabi¬ 
lity of this new profile was examined and its critical Reynolds 
number calculated. The new profile displayed a point of inflection 
near the wall. The new profile was found to be slightly less stable 
than the original Blasius profile. The minimum critical Reynolds 
number is reduced by about 52 in the example presented by Lessen 
and Gangwani. This amount of change is probably not very signifi¬ 
cant experimentally. 

Merkle, Tzou and Kubota [6] assumed that roughness enhances mo¬ 
mentum and heat transfer near the surface and modeled this using a 
turbulent eddy viscosity in the vicinity of the roughness, in a 
region they called the roughness sublayer. They concluded that the 
roughness sublayer alter? the mean velocity and temperature profiles 
within the laminar portion of the boundary layer. They obtained an 
increase in shape factor primarily through increase in displacement 
thickness. Stability calculations for this distorted mean profile 
indicated reduced stability. An attempt by the authors to validate 
their results by comparison with Achenbach's experimental results 
for a circular cylinder [7] resulted only in qualitative agreement. 

Experimentally, there is as yet no demonstration that for dis¬ 
tributed roughness, a profile change is responsible for enhanced 
instability. The present work intends giving a definitive experi- 





mental look at this question. Accordingly, the mean velocity pro¬ 
files, disturbance energy distributions and their spectral develop¬ 
ment are studied in order to help understand the mechanism by which 
distributed roughness influences transition. The investigation is 
conducted on a flat plate model with two roughness sizes in addition 


to the reference case of a smooth surface. 


CHAPTER II 


TEST FACILITIES 

2.1 Low Spe ed Mind Tunnel 

2.1.1 Physic a l Desc ri pt ion 

The low speed tunnel is of the open circuit type, with a con¬ 
traction ratio of 10.4:1 (Fig. 3). There are two screens ahead of 
the contraction, each with an openness of ft5%. The tunnel is ope¬ 
rated by a 53 IIP DC axial fan. The test section is 7.5' long and 
its cross-section is 28" x 28" (Fig. 4). Except for the bottom wall, 
the rest of the test section walls are made of transparent plastic 
(plexiglas). 

2.1.2 CjtJ i brat ion 

The inside of the tunnel was cleaned and sealed to achieve as 
low a froe-stream turbulence as possible. At the plate location, 
values of free-stream turbulence intensity below 0.03% at a test 
speed of f> m/s or were obtained as shown in Figs. 5 to 7. The hot¬ 
wire readings are affected by the velocity components which are per¬ 
pendicular to the wire; thus rite measured turbulence intensity is: 



2 • 2 Elat P late Model 
2.2.1 The Plate 

The model is a flat smooth plate, 41.5" long, 25" wide and 
i/8" thick (Figs. 8-10). The leading edge has a 1/32" radius fol- 
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lowed by a smooth taper which blends into the flat plate surface 
1.2" behind the leading edge. It is the same leading edge contour 
used by Strazisar et al. [8, 9] in their experimental study of the 
stability of heated laminar boundary layers in water. The surface 
of the plate was polished to a mirror finish using aluminum polish 
compound Met-All. No. 1187. 

2.2.2 The Mounting 

The plate is supported by two swept front supports which are at 
a 70° angle to the plate, and two rear supports perpendicular to the 
plate. The front supports have a symmetric airfoil cross-section 
with 2 3/4" chord, 3/4" thickness and the cross-section center at 
the plate junction is 16.5 cm downstream from the leading edge so 
as to prevent separation from the supports and generally minimize 
creation of additional disturbances in the tunnel. The rear sup¬ 
ports are 3/4" diameter cylinders. Their center is 80 cm downstream 
of the leading edge. They are not as important as the front sup¬ 
ports because their location is behind the test region. 

The choice of plate location was based on the wind tunnel ca¬ 
libration results. Since the longitudinal variation of free-stream 
turbulence intensity is constant, the leading edge of the plate was 
located near the door of the test section for convenience in making 
it easier to clean the plate surface. It was positioned at a 17.5" 
height above the floor and symmetrically in the horizontal plane. 

The natural vibrations of the plate on its supports were meas¬ 
ured at a tunnel speed of 6 m/sec. Fig. 11 shows that the highest 
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amp 1itudo appears near the leading edge and is 0.04 mm and occurs at 
a frequency of 8 Hz which is below the theoretical critical frequency 
hand for TolImien-Schlichting waves. Its magnitude is much below 
the roughness heights tested. The amplitudes at the middle and end 
of the plate are about half those near the leading edge. A flap at 
a 1° angle into the test-stream at the end of the plate was used to 
stabilize the stagnation point near the leading edge, in order to 
prevent leading edge separation. 

2.2.1 Roughness 

Two sizes of sandpaper made by the Norton Company were used. 

They were applied to the plate by using RTV adhesive. The finer 

sandpaper (80 grit) has a specified nominal average grain size of 

266 microns (0.0105 inch). Thus Re, at X = 30 cm at U =6 m/sec is 

k e 

10.6. A filler was applied starting 1 cm upstream of the leading 
edge of the sandpaper and blending smoothly into the sandpaper 
leading edge. 

The coarser sandpaper (24 grit) had a nominal average grain 

size of 1015 microns (0.0408 inch) so that Re, at X = 30 cm and 

k 

1) = 6 m/sec is 155. In this case, the filler began 1 inch ahead of 

the leading edge of the sandpaper. Because of the different thick¬ 
nesses of the two sandpapers, the ramp angle of the filler was 
approximately the same in both cases. 

Both sandpapers covered the whole plate surface downstream of 
their leading edges. For the small roughness, Re^ was calculated 


using the formula: 











Re k = 


= 0.332 (—) k X k 


( 2 . 2 ) 


based on the fact that for small heights above the plate 
vX -1/2 

—— = 0.332 U (——) iust as at the surface [1]. Re could also 

dy e U k 

R 

be obtained by calculating the n corresponding to k, n = kw —, and 

then the velocity u, could be taken from the Blasius curve. This 
k 

method was applied for the coarser sandpaper where the grains are 
not small. 

2.3 Measuring Equipment 

2.3.1 Pitot-Static Tube 

The pitot-static tube used was of the Prandtl type made by 
Sensor United, model PAC-12KL (Fig. 12). It is inserted from the 
top of the test section 30 cm downstream from the leading edge and 
displaced laterally from the hot-wire probe by 5 cm to minimize any 
interference effect. It is insensitive to yaw or pitch angles up 
to 30° and its quoted error is less than 2%. It has 1/8" outside 
diameter and is positioned 5 cm above the plate. Thus, the distance 
from the surface is more than 12 boundary layer thicknesses. Since 
the probe is also more than 10 probe diameters above the plate, the 
effect of the plate on the readings is negligible. 

The Reynolds Number based on the impact hole diameter (0.5 mm) 
for 6 m/sec speed is approximately 100. In such a case when it is 


larger than 30, viscous effects in the internal flow can be neglect¬ 
ed. The viscous effects in the external flow are also negligible 
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when the Reynolds number based on the probe radius is larger than 
100, which is the case here. 

For T.I. around 0.03% and dynamic pressures of about 0.1 in. 

1 2 

H.jO.Lhe Fluctuation velocity component dynamic pressure, — p(u') 

-9 

corresponds approximately to 9 x 10 in.l^O. This is negligible 
relative to 0.1 in.H^O. Thus, the pitot tube reading is insensitive 
to small turbulence. 

The time constant for the pitot probe is approximately five 
seconds. 

The probe is connected to a Meriam manoaeter which it ebl~ 
to measure 0.00 1 in.il ? 0 (Fig. 15). 

2.1.2 Hot Wire 

The shape of the hot-wire probe was designed appropriately for 
boundary layer measurements (Fig. 13). It consists of 1/4" stainless 
steel tubing in which electrical wires run through, that are con¬ 
nected to the copper prongs. The prongs are insulated one from the 
other by 5 minute epoxy. The tips of the prongs are connected by 
Wollaston wire (90 7. platinum, 10% rhodium) 8 mm long. 

The length of the control portion which is the sensor varies 
depending on the etching procedure during its preparation and is 
approximately 0.8 mm long. Beneath the joining of the prongs there 
is a protruding support made of 5 minute epoxv. From cathetometer 
measurements supported by mean velocity profile results, the wire 
is 0.2 mm above the support bottom. 

The principle of the hot-wire and the factors affecting its 
operation are described briefly as follows: 
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The wire is heated by a constant current. 1t is cooled bv the 
flow causing the temperature to drop and,consequently, the electric 
resistance to diminish. It is cooled primarily by forced convection 
and heat conduction. Total amount of heat transferred depends on: 

1) The flow velocity 

2) The difference in temperature between the wire and 
the fluid 

3) The physical properties of the fluid 

4) The dimensions and physical properties of the wire 
The last three are known and thus the flow velocity can be 

measured. 

The hot-wire anemometer is still a prime instrument for tur¬ 
bulence measurements. Its principal advantage is the separation of 
the mean velocity from the fluctuating components and the ability Co 
measure them separately at the same time. The other main advantage 
is in its being thin. The Reynolds number based on wire diameter 
is 1.34 for a speed of 6 m/sec. This is below 3, which is the limit 
for smooth flow around the wire without separation even with heating, 
because Reynolds number decreases with temperature | JO] (fig. 14). The 
thermal inertia of the thin wire is very small and thus the response 
for frequencies of interest in this experiment (up to 5 KHz) is 
practically instantaneous. It catches less dirt than a thick sensor 
and altogether it is sufficiently strong and rigid to sustain 
vibrations and stresses caused by turbulence. The measured velocity 
is in good agreement with the flow velocity as measured using the 
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pitot-static tube for turbulence intensity smaller than 10%. 

The joining of the prongs is 32,000 wire diameters behind the 
wire and thus there is no real influence of the prongs on the flow 
around the sensor portion of the wire. 

There are two problems witli the hot wire. The first is 
drifting and the second is contamination. Both can be treated prop¬ 
erly as will be explained in the next chapter. More about hot wires 
can be found in ref. {11 1. 

The hot-wire anemometer set is a Shapiro and Edwards model 50 
(Kig. 15). It includes an amplifier, a resistance bridge, a poten¬ 
tiometer, a mean square output meter and a square wave generator. 

The potentiometer allows measurement of wire currents up to 10 ^ amp. 
and so enables accurate velocity measurements. A compensation cir¬ 
cuit is available but is needed only for frequencies above 5KHz and 
was not utilized in this experiment. 

The preparation of the hot-wire sensor is as follows: First, 
the wire is soldered to the prongs and then the probe is attached to 
a positioner. The wire is brought into contact with a small drop of 
707.', concentrated nitric acid ( p igs. 16-17). An electrolyzing cur¬ 
rent is supplied bv a 1 1/2 volt Eveready battery for about one 

minute. Tiie length of the uncovered (etched) center portion of the 
wire is a function of the drop size, and is of the order of 0.8 mm. 
This portion which is the sensor has a diameter of 0.00015 in. 

The wire is then checked using a microscope which magnifies 
60 times. It is important to tiave a straight and clean wire. 
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The soldering of the wire to the prongs with a soft solder had 

to be done quickly in order not to create excessive heat in the wire 

and thus avoid any elongation. If overheated, the uncovered portion 
of the wire could not sustain as much stress as the rest of the wire, 

and at the end of the etching process would become deformed. 

The other factor which affects wire shape is the surface tension 
of the drop. To avoid any adverse effect, the wire has to be taken 
out of the drop very carefully. 

2.3.3 Traverse Mechanism 

The traverse mechanism (Fig. 13) was designed to give accurate 
vertical readings of the probe position. It has a metric Mitutoyo 
micrometer which can be moved longitudinally with the help of a 
motor. There is an error of ^0.2 mm in longitudinal positioning. 

The micrometer enables vertical precision to '0.01 mm and so this is 
the accuracy of the hot-wire probe position because it is attached 
to tlie micrometer. The height can be read on the micrometer scale 
and readjusted to zero each time the hot wire is brought to the 
plate surface. The precision of lateral positioning is not as impor¬ 
tant as the vertical one and there is an estimated error in lateral 
position of at most 0.5 mm. 

The traverse mechanism is outside the top of the test section 
and the hot-wire probe goes through and moves in the slot made in 
the top wall. This location of the mechanism was chosen because it 
allows the probe to have a simple shape, has a shorter shaft than 
a probe from the side wall and does not vibrate at the test speed. 
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(t also creates less disturbance than a probe from the side wall. 

The slot was placed above the middle of the plate where the 
turbulence intensity is smallest. 

2.1.4 Spectrum Analyzer 

A Unigon Model 4512 spectrum analyzer was used to measure 
spectra. It has an averaging capability for a chosen time interval. 
The bandwidth and the amplitude scale could be adjusted. A Polaroid 
camera attached to the screen enabled pictures of the spectrum 


to be taken. 





CHAPTER III 


EXPERIMENTAL PROCEDURE 


3.1 Ca 1 jLb ration 

Calibration of the hot wire was performed with every run 
because the cold wire resistance was changing due to heat, dust and 
natural deterioration. To eliminate drifting, it was necessary to 
allow two hours for the system to warm up before using the hot wire. 
The contamination problem was solved by changing the wire every few 
weeks. 

Calibration was always performed with the wire at 3 cm above 

the plate, and 30 cm downstream of the leading edge. The cold wire 

resistance R was measured. The overheat ratio of 1.3 was ad- 
o 

justed by setting the wire resistance R to a value of 1.3 R . 

w o 

The wind tunnel speed was increased gradually and each time the 

wire current and the manometer readings were taken. In such a 

manner, the calibration curve for the velocity as a function of the 

wire current was obtained and plotted (Fig. 19). It was almost a 

straight line as it should be according to King's law [111: 

L 2 = a + b/u (3.1) 

w 

After the calibration, the pitot tube was pulled out. 

3.2 Smooth Plate 

3.2.1 Adjustment of Zero Pressure Gradient 

The outer flow velocity was measured along the plate using the 
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the critical stability region; that is where, according to stability 
theory, the disturbances are amplified. Since tile concern was to 
achieve a constant pressure in this region, the rear supports were 
shortened by I cm, and so the plate was tilted by a 1° angle. The 
angle of attack could he control Led further by elevating the rear of 
the plate using the screws in the rear supports but this was not 
necossa ry. 

As can be seen from Fig. 22, this adjustment yielded a constant 
pressure region between X = 30 cm and X = 60 cm. The mean velocity 
profiles in this region show good agreement with the theoretical 
Kins ins curve (Fig. 23). In this region, the Reynolds number based 
on displacement thickness is 390 to 83 r ' which is in the critical 


stability range for amplification of Tollmien-Schlichting waves. 
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Thus, the region between X = 30 cm and X = 60 cm was chosen as the 
test region. 

3.2.2 Measurement of Mean and Fluctuating Information 

The mean and the fluctuating information was taken at the same 
time. For the smooth plate, the probe was lowered until the support 
touched the plate. This was repeated several times and the micro¬ 
meter readings were taken. They were always within r0.03 mm of a 
mean value and thus could be considered accurate enough. At this 
position, the wire was approximately 0.2 mm from the surface. This 
was the closest position from the plate and there the first readings 
were taken. Measurements were repeated at vertical intervals of 
0.3 mm until there was no current change, indicating the outer flow 
regime. 


At each point, the current, the series resistance and the 
average wire voltage V(e ')^ were measured. This enabled calculation 
and extraction of the information of interest. The mean velocity 
was taken from the calibration curve (Fig. 19) for the corresponding 
current. The fluctuation rms velocity component was calculated as 
follows [111: 


(u 1 ) 


2 = 


U 


e') ‘ 


(3.3) 


R I aEZ 
w w 

where S =-is a hot-wire sensitivity 

2u 


R = 1.3 R - the wire resistance at overheat 
w o 

I - 
w 


the wire current 




a 


= (K -R )/R = 0.3 

w o o 

!■; = (1-. )/(l+2'u) where r = R /(R + R^) and 

H - the series resistance 

1/(1+ /u /fn) where H is the corresponding velocity 
o o 

for a zero current and extrapolated from the cali¬ 
bration curve 
,, - the mean velocitv 


The volt.iye '/( e') “ was measured by taking approximately 50 


readings at 5 second intervals. 


By calculating 


lurbance energy 


;y Vo *') 2 £ 


n/=f 

V(u') = u , it was easy to get the dis- 

rms 


and the free-stream turbulence intensity 


V<V) 2 /u 


or the local intensity 


V(u’) 2 /u. 


1. 1 Ro ugli Pl ate 

1. 1.1 Ad ditional Adjustments 

The finer sandpaper was at first positioned to begin at 3 cm 
downstream of the leading edge. The Re^ and k/A* at this location 
are 33 and 0.56 respectively. The boundary layer was apparently 
tripped by the leading edge of the sandpaper and the flow at the 
l irst measuring station,X = 30 cm,was already turbulent. The paper 
was then cut back to 6.5 cm downstream of the leading edge of the 
plate. The corresponding Re^ and k/j’* were 23 and 0.38 respecti¬ 
vely. Here the flow was laminar at the first station and for this 
location of the leading edge of the sandpaper, the measurements 
were taken. 

The coarser sandpaper began initially at 6.5 cm from the lead- 
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ing edge of the plate but again the boundary layer was tripped. 

The corresponding Re^ and k/6* were 314 and 1.48 respectively. 

The flow was turbulent even at X = 20 cm. When the paper was cut 
back to 18.3 cm, the flow at the first station was laminar. At 
X = 18.3 cm, Re^ = 192 and k/S* = 0.88. This was the location of 
the leading edge of the coarser sandpaper for the measurements 
presented herein. 

3.3.2 Measurements Taken with Each of the Tw o Sandpapers 

The measurements and the procedure for the finer sandpaper was 
similar to those for the smooth plate. For measurements with 
coarser sandpaper, precautions had to be taken to prevent wire 
damage. 

A glass plate of 1.4 mm thickness was placed on the surface 
under the probe and the probe was lowered until the support touched 
the glass plate. The micrometer reading was taken. The probe was 
raised and the glass plate removed. The probe was lowered again by 
1.4 mm beneath the reading. At all stations, it was possible to 
lower the probe an additional 0.5 mm without damaging the wire. The 
probe location was illuminated, and it was always possible to see a 
gap between the wire and the grit. This was checked many times. This 
is possible because of the non-uniform distribution of the grains 
in the sandpaper matrix and also because of some grains which were 
positioned above the average. If the plate rested on those higher 
grains, then the wire could be lowered beneath their height. We 
could assume that there were enough grains of such size to allow 
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tii is to occur. Tit us, the first reading for the roughest sandpaper 
was taken with the wire at 0.3 mm below the highest grains in the 
near vicinity where the measurements were taken. 

In addition to the measurements as taken earlier, special 
data were taken for the coarser sandpaper. Here the information 
was taken also at X = 20 cm. The hot-wire probe was positioned at 
the peak energy location. The amplitudes were measured from Pola¬ 
roid pictures of the analyzer output and the amplification rates 
were calculated from: 


a. 

l 


1 dA 
A dX 


(3.4) 

The amplitude gradient was measured from the amplitude curves. 

Data were also taken using the filters on the hot-wire anemo¬ 
meter panel. At the peak energy location for each station, the 
energy was measured for different frequency bands. The filters were 


1-10, 10-100, 100-5000 Hz. 






CHAPTER IV 


RESULTS AND DISCUSSION 


4.1 S mooth P la te 

To provide a reference for the later measurements with surface 
roughness, a set of mean and fluctuating measurements was taken 
on the smooth aluminum plate. The mean velocity profiles at X = 30, 
40, 10 and 60 cm downstream of the leading edge are shown in 
ligs. 24 (a-d). Two data sets taken at different times are shown. 

The Blasius curve taken from Sehlichting [12] is plotted for com¬ 
parison on each of the figures. The points agree with the Blasius 
curve within the experimental error. The profiles corroborate 
the earlier indication from surface pressure measurements that the 
region 30 • x 60 cm is a constant pressure region. 

One set of test results for the fluctuation energy nondimension- 

2 

a Iized with respect to U~ versus dimensional boundary layer coordi¬ 
nate v and nnndimonsional n is presented in Tigs. 25a and 25b 
respective1v. 

The free-stream turbulence and the wind tunnel fan rotation 
affect the magnitude of the disturbance energy, but the shape and 
order of magnitude of the curves remains consistent. Even though 
the data points are not always on the curve lines, the deviation is 
within the experimental error and the curve lines are definitely 
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rharaeteristic of the phenomena taking place according to other 
test results taken in the present investigation. The shape is si¬ 
milar to that obtained in other experimental studies for a smooth 
plate, for example that of Strazisar et al. [8,9]. 

The integrals under the curves and the peak values are increas¬ 
ing consistently and gradually in the downstream direction. The 
energy content is distributed over the whole boundary layer but is 
concentrated mainly in the middle. It is almost symmetrically dis¬ 
tributed. The peak moves away from the wall but remains at 

n 2.3 or y /<*> ~ 0.4 3 and y /6 *1.37. The peak energy is 
max max 

increasing approximately linearly with distance downstream. 

The local disturbance intensity which is with respect to the lo- 
caI velocity is plotted in Tig. 26. It is seen to decrease near the 
wall witli downstream distance. The energy integrals and peak values 
are increasing but very little. Both the energy and the local 
intensity are increasing in the vicinity of the edge of the boundary 
layer, thus indicating that the response to free-stream disturbances 
is not confined to tiie interior of the boundary layer. The turbu¬ 
lence intensity at the outer edge is seen to increase about 0.02% 
with every 10 cm downstream from the leading edge. 

The reasons for showing both energy and local turbulence inten¬ 
sity will become apparent later in the cases with roughness. 

h.2 Distributed Roughness of Re^ (X = 30 cm) = 10,6 

Similar measurements were taken for this roughness as for the 
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smooth surface. The mean velocity profiles at X = 30, 40, 50 and 
60 cm downstream of the leading edge are shown in Figs. 27 (a-d). 
dust as for the smooth case, a number of data sets were taken and 
the Biasius curve is plotted for comparison on each of the figures. 
The zero reference height, as for the smooth wall case, was assumed 
to be 0.2 mm beneath the hot-wire when the probe support touched 
the surface. With this choice of n = 0 (y = 0), the Blasius curve 
fits the data points within the experimental error. 

One set of test results for the energy is presented in 
Figs. 28 a,b. The shape and order of magnitude of the curves for this 
run remains consistent for the results of other tests taken in this 
investigation but not shown. The energy has the same shape and the 
same gradual increase as for the smooth case. It is again mostly 
concentrated at the center of the boundary layer and the peak 
remains at n ' 2.3 as for the smooth surface. 

Wo can see already slightly more curvature near the wall than 
in the smooth wall case (Fig. 25b) as the flow proceeds in the 
downstream direction. The energy at the outer edge of the boundary 
layer is also slightly larger than for the smooth plate. Within the 
boundary layer the energy levels are slightly smaller than for the 
smooth surface perhaps because the speed was lower by 0.2 m/sec and 
possibly because the flow was less turbulent initially, thus affecting 
the initial energy levels. Although the energy peak at X =30 cm is 
smaller than for the smooth wall case, the energy growth - the 
ratio of energy at the last station to that at the first station - 






is 107. larger than before and that, perhaps, is an indication of a 
slight roughness effect. 

A more pronounced indication of the roughness effect is shown 
by the local turbulence intensity distribution presented in Fig. 29. 
Even though the energy was smaller, the intensities were larger than 
those for the smooth wall (Fig. 26) in the vicinity of the wall. 

The intensities of the peaks close to the wall are increasing with 
downstream distance in contrast to the results for the smooth 
surface. 

In summary, this small roughness (beginning from Re^ = 22.8 
and k/(S* = 0.38 at X = 6.5 cm) has only a very small effect. It 
doesn't perturb the flow to the nonlinear regime and both the mean 
velocity and the disturbance energy profiles resemble those of 
laminar flow on a smooth plate. Although the overall effects are 
small, we could detect changes in the local turbulence intensity 
distributions, particularly in the vicinity of the wall. 

4. 1 Distributed Roughness of Re^ (X = 30 cm) = 155 

4.3.1 Mean an d Fluctua tion Meas urements 

A set of mean velocity profiles is shown in Figs. 30 (a - d). 

In this case the roughness is large enough so that the effective 
zero velocity location for the boundary layer could not be deter¬ 
mined a priori. The data point closest to the wall was 0.3 mm 
below the Lop of the largest grains. This point was chosen as the 
reference location for the rest of the data taken. The nominal 
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gri tsize is known and thus the distance of the reference location 
point from the nominal bottom of the grains is also known. By 
fitting a curve through the data points, the extrapolated zero velo¬ 
city point can be located. For the firs't two stations (X = 30 cm, 

X = 40 cm), a Blasius curve could reasonably be put through the 
data points and the extrapolated zero velocity locations were at 
approximately the same distance from the nominal bottom of the grit 
(0.2 mm). An attempt to fit the data points with a similarity pro¬ 
file for 8 = -0.1 was unsuccessful. These two measured profiles are 
much closer to the Blasius curve and compare very well with the 
measured data on the smooth plate once an adjustment is made for 
the extrapolated zero. 

Curves through the mean velocity data points for the two down¬ 
stream stations (X = 50 cm, X = 00 cm) are shown in Figs. 30c and 
30d. The extrapolated zero velocity location is chosen to be the 
same as in Figs. 30a and 30b. Blasius curves are drawn for compa¬ 
rison. At X = 50 cm the mean velocity profile has changed signifi¬ 
cantly from that at X = 40 cm. At X = 60 cm there is the beginning 
of the emergence of a turbulent profile shape. 

if there are inflection points near the wall in Figs. 30a and 
30b, it is practically impossible to detect them considering the 
experimental error. Tt is also very difficult to Jo a detailed 
search in the vicinity of the grit because of prospective wire 
damage and probe interference effects. 
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Fluctuation data for two separate tests taken a week apart are 
shown in Figs. 31a, b. The first test results are indicated by I 
and the second by II. The first test corresponds to the mean velo¬ 
city results shown earlier (Fig. 30).* The n = 0 locations used in 
these figures for each station were taken from the extrapolated zero 
velocity points of the mean velocity profiles. The shape of the 
energy distributions is not smooth but somewhat oscillatory. Both 
tests show more than one local peak and the number increases with 
the distance downstream. The first two stations display more or less 
their previous shape. Their peaks are in the vicinity of n * 2.3 
and only the second test deviates from this slightlv at the second 
station. The downstream profiles here are much fuller than for the 
smooth plate or the finer roughness and there is significant 
increase in energy content near the wall. 

For the finer roughness the only clear indication of a rough¬ 
ness effect was through the changes in the turbulence intensity dis¬ 
tributions. Here we see a large effect on the energy distributions 
as well. The peaks are moving closer to the wall and most of the 
energy content is concentrated in the n < 4 region for the first 
three stations. Farther downstream there is also significant 
increase in energy content at the outer edge of the boundary layer. 


* Mean velocity profiles for test II (4/21/80) were measured at 
X - 30 cm and X = 40 cm. they are verv close to the profiles of 
Figures 30a and 30b. The profiles at X = 50 cm and X = 60 cm differ 
slightlv from those of Figures 30c and 30d in a manner consistent 
with the different progressions to transition on those two days. 
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In between X = 40 cm and X ~ 50 cm whore a velocity profile change 
takes place, the energy increases by an order of magnitude. Figs. 
32a,b present the local turbulence intensity for these two tests 
respectively, again using the extrapolated zero from Test 1. The 
local intensity rises to more than 15% between X=A0cm and X = 50 cm 
emphasizing that the disturbances are in the non-linear regime. 

In Test I, there is not much overall energy change between the 
last two stations but there is a noticeable buildup in the outer 
portion of the boundary layer. Accompanying this is a definite 
change in the mean velocity profile (Figs. 30c,d) and an increase 
in the boundary layer thickness even though it is not yet a fully 
turbulent flow. 

Test II shows a flow that seems less disturbed initially. Less 
of an increase in energy occurs between X = AO cm and X = 50 cm. 

As a consequence, there is still an increase in energy between the 
the last two stations. Only the last station indicates a signifi¬ 
cant energy increase in the outer portion of the boundary layer. 

The mean velocity profiles which are not presented here, again 
indicate a Blasius profile for the first two stations and a very 
slight change relative to Test f at the third station. The inten¬ 
sity peak at the last station reaches the same magnitude as at the 
last two stations of Test I. 

The intensity shape is bumpy. However, the main peaks are approach¬ 
ing the wall t consistent with the observation for the Re^ = 10.6 case. 
The peaks are within the region n •- 0.5. At the measurement point 
closest to the wall, the total amplitude ratio relative to free- 











stream turbulence level (assumed to be 0.02%) for both experiments 

7.3 

is already about 1500 which is c . The boundary layer at this 
point is however not yet fully turbulent. Maybe the ratio is larger 
closer to the wall hut it was impossible to measure any closer. 

The consistency of the peak intensity development shown by 
Figs. 32a, b with that of the finer roughness indicates that the 
extrapolated zero point is reasonably chosen. The development of 
the energy and the intensity indicate that the fully turbulent 
regime is not yet reached. 

An attempt to isolate a roughness effect for Test II by 
subtracting Fig. 26 (the smooth plate intensity) from Fig. 32b (the 
rough plate intensity) is presented in Fig. 33. The curves are 
essentially the same as in Fig. 32b indicating that the overwhelming 
contributions to the data in Fig. 32b are due to the roughness. 

Fig. 34 shows the variations of the peak local intensity and 
the square root of the peak nondimcnsional energy with Re^. These 
results are consistent with what was said earlier: (a) the intensity 
is more sensitive than the energy and thus a better indication of 
the degree to which the flow is disturbed; (b) the flow in Test I 
Is more disturbed. 

In summary, the intensity magnitudes for both experiments are 
above 15% at the third station (X = 50 cm) and thus are clearly non¬ 
linear. The intensity for the second station is 4.5 - 8% and thus 
significant too and probably non-linear. The main conclusion is 
that the disturbances grow to non-linear amplitudes before detect¬ 
able velocity profiles are observed. 








4.3.2 The Spe c trum 

In addition to mean and fluctuation measurements, one set 
of tests was made using the llnigon spectrum analyzer. Figs. 33 - 39 
show the spectra taken at X = 20, 30, 40, 30 and 60 cm with the 
coarsest paper and Li ^ = 6 m/sec. For X = 20 cm to X = 40 cm 
(Figs. 35, 36, 37) gain settings were inadequate to get data above 
20 Hz from those pictures. It was possible to read tile amplitude 
above 20 Hz for the last two stations but the accuracy is doubt¬ 
ful, and since the amplitude slope is very sensitive to any small 
amplitude change, the calculated v values are doubtful too. 
Nevertheless it is possible to extract some important information 
from this test. 

Amplitude distributions for frequencies up to 17.2 Hz are shown 
in Fig. 40. They have a shape resembling what is obtained in linear 
stability studies. Fig. 41 presents the calculated amplification 
rates For 17.2 Hz amplification occurs at Re . - 330. As the 

i 6 * 

Re increases, successively lower frequencies become amplified but 
their maximum amplification rate is diminishing. The highest ampli¬ 
fication rate for 17.2 Hz is located at Re , = 650. Lower frequen- 

<5* 

cies show peak rt . at successively larger values of Re . These 
1 (S * 

peaks occur in the non-linear region of the flow as identified from 

energy and intensity distributions. Fig. 42 shows again that as 

Re ^ increvises, the lower irequeueies become important. A plot of 

“j vrsus Re. A (Fig. 43> indicates lhat higher frequencies (above 
max 
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17.2 Hz) may be important below Re = 650. Note that the Tollmien- 

o* 

Schlichting hand at Re . = 650 is between 45 - 85 Hz. In fact if 

o * 

we extrapolate the trends of Fig. 41 to higher frequencies, it would 
suggest the higher frequencies become more important close to the 
leading edge. 

Ampl if icat ion rates -r*. for the frequencies up to 130 Hz at 
X = 50 cm are shown in Fig. 44. Similar behavior of a. is found at 

i 

X = 60 cm. Amplification is shown at all frequencies. There does 
not seem to be anything special about the Tollmien-Schlichting 
hand. Again, at X = 50 cm, the disturbances are non-linear. 

The experimental "neutral curve" for frequencies up to 17.2 Hz 
is below the neutral curve for the linear stability of non-parallel 
flat plate flow as calculated hy Saric fl3] (Fig. 45). It is also 
below the neutral curve for the linear stability of non-parallel 
I’a 1 kner-Skan flow corresponding to 6 = -0.1. This corroborates the 
earlier indication that disturbance growth was governed by non¬ 
linear growth effects rather than by linear instability of an in¬ 
flectional mean velocity profile which hasn't been detected in any 
case. 

From the trend of the curve we can definitely say that there is 

earlier amplification and very likely earlier transition than for 

the smooth surface. While the linear theory for a smooth plate 

predicts :l . for Re ^ above 684, here it is larger much below this 
ma x 

number. The order of magnitude of (a. 6* ) is also totally 

1 r r ma x 


dif ferent . 
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-a. 6* 

The maximum value of ( — yt ) for non-parallel smooth flat 

66 -6 

plate flow is approximately 9.55 x 10 and occurs at Re^ = 795 [9]. 
For the same Rethe present data for the rough case shows it equal 
2.67 x 10 ^ and the value of is even higher for lower Re^. This 
means faster transition. Considering the state of the boundary 


layer at X = 60 cm (Re^ = 835), transition might be completed at 

Re ' 350,000 which is very low compared to the value of 3 x 10^ 

X 

that might be expected on a smooth plate for the low free-stream 
turbulence intensity of the tunnel. 

Frequencies corresponding to roughness wavelengths are of order 


— = (~~) w • ln the present experiment such frequencies would be in 
the range of 1.5 to j KHz. No growth was observed in these frequency 
bands. 


4.3.3 Fi1 to r ing 

Another experiment was performed using the built in filters on 
the hot-wire anemometer set at the peak energy location at X = 30, 
40, 50 and 60 cm downstream from the leading edge. The filters were 
l — 10, 10-100, 100-5000 Hz. Fig. 46 shows the energy results for 


these frequency levels. 

There is a shift in energy to higher frequencies. Most of the 
energy in the non-linear regime is concentrated in the 100-5000 Hz 
band. This didn't show up from the spectrum because even though 
the amplitude in the high frequency range is small, its integral is 
significant and very sensitive to a small amplitude increase, since 
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the band width of the 100-3000 Hz band is so large. Major high 
frequency content is probably in the 100-200 Hz band. The mechanism 
is driven by frequencies up to about 140 Hz as indicated by the 
spectra (Figs. 3d, 39), but the energy is dissipated into higher 
f requencies. 

4.4 Lea ding E dge of the Dis tr ibutcd Roughness 

As described in the preceding chapter for both roughness sizes, 
the leading edge of the sandpaper had to be cut back in order that 
the boundary layer not be tripped by the leading edge of the sand¬ 
paper. Such tripping did occur before cutting back the sandpaper, 
despite the use of a filler to avoid a step at the leading edge of 
the sandpaper. 

The finer roughness sandpaper was cut from X = 3 cm with the 
corresponding Re = 33 and k/<S* = 0.56, to X = 6.5 cm downstream from 

K 

the leading edge and the corresponding Re^ = 22.8 and k/6* = 0.38. 
Before the cutting hack the sandpaper, the flow was turbulent over 
the whole test region and after cutting it back it was laminar. 

The coarser roughness paper was cut back from X = 6 cm with the 
corresponding Ke^ = 314 and k/<S* = 1.48,to X = 18.3 cm with corres¬ 
ponding Re = 198 and k/fi* - 0.88. Before cutting the sandpaper back, 
transition was detected already at X = 20 cm and turbulent flow oc¬ 
curred in the whole test region. After cutting back the sandpaper, 
laminar mean profiles were obtained at least to X = 40 cm (Fig. 30). 








Eiven though Re^ was equal to 198 at the leading edge of the 
coarser sandpaper and only 33 at the leading edge of the finer sand¬ 
paper and the filler ramp angle was approximately the same in both 
cases, the latter tripped the boundary layer while the former did 
not. Thus Ke is not the only parameter affecting the flow. The 

K 

same can probably be said about k/o*. Thus, the response of the 
flow to the leading edge of a rougness region cannot be simply 
described in terms of these two parameters. The filler strip 
geometry and the smoothness achieved are important factors in 
influencing the flow. 


I 









CHAPTER V 


SUMMARY AND CONCLUSIONS 

A preliminary experimental study has been completed of the dis¬ 
turbances in a laminar boundary layer due to distributed surface 
roughness. Tests were conducted on a smooth plate and with two dif¬ 
ferent grades of sandpaper covering. 

Over the portion of the smooth plate displaying zero pressure 

/ 

gradient (30 cm < x < 60 cm), the mean flow velocity profiles agree 

with the Blasius profile within the experimental error. The u' 

rms 

distributions display a peak at n * 2.3 that increases slightly with 

distance downstream reaching a level of 0jOO 7 U g at X = 60 cm. 

For distributed roughness having Re^ ® 10, the results are quite 
comparable to those for the smooth plate. The mean flow velocity 
profiles agree with the Blasius curve within the experimental error 
and the u^.^ distributions are closely comparable in level and 
distribution to those for the smooth plate. There is however indi¬ 
cation of a slight increase in local intensity in the very near vici¬ 
nity of the wall. It is well to point out that although the dimen- 

2 

sionless turbulent energy (u' /U ) and the local turbulent inten- 

rms e 

sity ( u, rmg / u ) involve the same fluctuation amplitudes, the local 
turbulent intensity Is a much more sensitive indicator of effects of 
roughness. 
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For distributed roughness with Re is 150, the picture is dif- 

K. 

ferent. At X = 30 cm and at X = 40 cm, the mean velocity profile is 
best fit by a Blasius profile whose effective origin is somewhere 
within the roughness height. Beyond X = 40 cm the mean profile 
seems to progress toward a turbulent profile. Up to X = 40 cm, the 
u' still peaks at n ; 2.3 but the local intensity is now of the 
order of 5%. Beyond X = 40 cm there is a large increase in the u' 
levels and an increase in the number of peaks. Very large local 
intensities are obtained in the near vicinity of the wall (n < 0.5) 
but there is also noticeable buildup over the entire boundary layer 
to the outer edge. A study of disturbance spectra for this case 
shows largest amplitudes and amplifications at frequencies well 
below the Blasius neutral curve. The maximum amplification rates 
are about three times as large as those for Tollmien-Schlichting 
instability over a smooth plate. Transition occurs at 
RCj. '• 350.000 which is much below the value of 3 x 10^ or so that 
might be expected for a smooth surface at the turbulence level of 
the tunnel. 

To conclude, one observes that disturbances can grow to non¬ 
linear magnitudes before velocity profile changes are observed. 

The indicated importance of frequencies below those of the Blasius 
neutral curve was not expected and requires more careful examination. 
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Sdicm.it ic diagram of effect of distributed roughness on 
transition on a f]at plate. 



The mean velocity distortion due to 2-D roughness 
elemont. 









































Tlu* front viow of tlu* model. 
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Figure 14. Flow Regimes About Unheated Circular Cy¬ 
linders of High Length-to-Dlame ter Ratio 
Versus Reynolds Number (Morkovin, 1964). 
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20. Pressure cot-f t ic ionc for tlio olate 
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f-'ig. 24c. Mean velocity profile for a smooth plate, X = 50 cm. 
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Fig. 27b. Mean velocity profile for a rough plate, Re (X = 30 cm) = 10.6; X = 40 cm. 
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profile for a rough plate. Re. (X = 30 cm) = 153; X = 40 cm 




















































Tost I - Loc.il turbulence intensity profiles for a rough plate, 
Ke (X = TO cm) = 135. 
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II - Turbulence intensity profiles for a rough plate 
X = 10 cm) = 165. 
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li'.sL (I - boc.il turbulence intensity due to roughness only, 
Ko k (X = 30 cm) = 15 r >. 
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The liif;lu'st amplification rate versus X. 
Ue, (X = 30 cm) = 155. 
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Non-para 1 le 1 theoretic.il neutral-stability 
curve for ( =-0.1, Saric (1975) 

Non-parallel theoretical neutral-stability 
diagram for (• = 0, Saric (1975) 

Roughness results for (■ = 0, Re^ (X = 30 c 
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Neutral stability results. 





















